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In this paper a 4-component Dirac theory of ionization of hydrogen molecular ion in a super- 
intense laser field is presented. Simple analytic expressions for the spin specific as well as the total 
ionization currents emitted from the ground state of the ion are derived. The results are given for 
all polarization and finite propagation vectors of the field. They apply for the inner-shell ionization 
of analogous heavier molecular ions as well. The presence of molecular two-slit interference effect, 
first found in the non-relativistic case, and the spin- flip ionization current, and an asymmetry of the 
up- and down-spin currents similar to that predicted in the atomic case, are found also to hold for 
the present relativistic molecular ionic case. Finally, the possibility of controlling the dominant spin 
' currents by simply selecting the handedness of a circularly polarized incident laser field is pointed 
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PACS numbers: 33.80.Rv, 42.50.Hz, 32.80.Rm 



I. INTRODUCTION 

A non-perturbative relativistic analysis of intense-field ionization dynamics becomes necessary if: (a) the laser 
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i— | ■ intensity is so high that the ponderomotive energy Up = 4r ^ $ becomes comparable to or greater than the rest-mass 
energy m c 2 of the electron or U p > m c 2 , (b) the spin degrees of freedom, the magnetic field component and/or the 
♦p* retardation effect of the field become of interest, and (c) the initial bound electron is in a deep- lying shell of the target 
molecule or its ions, where the electron motion is already relativistic. ft should be also borne in mind that although 
initially a bound electron may be in the upper valence shell - where its motion can be non-relativistic ((3 — ^ << 1) 
-t— > \ - a super-intense laser pulse could, nevertheless, strip the target molecule of the outer shell electrons before the peak 
of the pulse arrives. Thus effectively the pulse at its highest intensities might in fact interact with the relativistically 
moving electrons of the inner-shells of the stripped ions, and ionize them further. Also, an initially non-relativistically 
moving outer-shell electron on virtual ionization at low velocities might be accelerated to relativistic velocities by 



the ultra-intense laser pulse. Furthermore, intense /lig/i-frequency free electron laser (FEL) radiations can interact 



directly and strongly with the relativistically moving bound electrons in the deep-lying shells of atoms and molecules. 

In view of the recent progress in ultra-short laser pulse techniques and high-frequency FEL technology, super-intense 
laser fields reaching intensities > 10 22 W/cm 2 , and FEL radiations of intensities > 10 16 W/cm 2 are now available or 
expected to be available in the laboratories in the near future. The purpose of this paper is to present for the first 
time a 4-component Dirac theory of ionization of the simplest of molecules, H^~, and of the analogous but heavier 
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K-shell molecular ions, by laser fields of arbitrary intensity and polarization and wavenumber. The present method of 
analysis within the strong-field KFR ansatz [3, 0] , is a generalization of the method developed earlier for the atomic 
case Q and can be further generalized to the inner-shell ionization of polyatomic molecules in a manner similar to 
the non-relativistic approach in[^, 



II. THEORETICAL FORMULATION 

The Dirac Equation of a diatomic molecular-ion interacting with an intense electromagnetic field can be written as: 

ih^ = [ca ■ ( Pop - * A(arf - « . r) ) - M _ M + ^ffo + ^ (1) 
at c r a rb R 

where, mo and eo are the rest mass and charge of electron, a and (3 are the standard Dirac 4x4 matrices, and 
Vi(t,r) = — eoct • A(u), with u = (tot — k ■ r), is the electron-laser interaction Hamiltonian and A(u) is the vector 
potential of the external radiation field: 

A(u>t — k ■ r) = Ao(ei cos (£/2) cos (tot — k ■ r) — e2 sin (£/2) sin (cot — k ■ r)) (2) 



We consider the general case of a laser field of arbitrary elliptic polarization (ellipticity parameter £[0,±7r/2]; ± 
stand for the left and right helicity) with the semi-major and the semi- minor unit polarization vectors e\ and e2, 



respectively. It should be noted that the Coulomb fields of the "heavy" nuclei a and b are treated as external fields 
acting on the electron at the position r. We have defined, r a = \r — R a \ and rt, = \r — Rb\; R a , Rb are the locations 
of the two nuclei a and b, respectively. Clearly, for the homonuclear diatomic case (CM. of the nuclei as coordinate 
center), R a = Rb = R/2 where R is the nuclear separation, and Z a = Zb = Z is the nuclear charge. 
We introduce the initial state partition of the total Dirac Hamiltonian 

H D (t)=H mol +Vi{u) (3) 

where, Vi(u) = —e$a ■ A(u), and 
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H mo i = [ca ■ p op H — - + m c (4) 

r a n R 

is the molecular Dirac Hamiltonian, with stationary eigenfunctions satisfying the Dirac equation: 

{E i -H mol )ip i {r) = Q (5) 

with positive and/or negative energy eigenvalues JSj. We choose the final state partition of the total Dirac Hamiltonian 
H~£>(t) (given by the right hand side of Eq. (JTJ) ) as: Hu(t) = Hy i{t) + Vf where, 

H Vo i(t) = (ca ■ (p op - —A(u)) + m a c 2 ) (6) 

is the Dirac- Volkov Hamiltonian of an electron in the plane wave field, Eq. and 

Vf - Vmoi{r; R) - — — + — - — . (7) 

'a 1 b 

For ease of working with Dirac equations it is useful to reduce the number of 4 natural constants appearing in them 
to only two by introducing the reduced constants: m — ' m B £ and e = f£ (which can be therefore restored easily in the 
end, if desired), and the following Feynman notations for 4- vector, 4-scalar product and the 4-gamma matrices: 
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(8) 

In this notation, for example, Eq. ([1]) takes the form: 

(< - e 4(x) - j V mo i{x; R) - m)V{x) = 0. (9) 

The total Dirac wavefunction can be expanded systematically, as in the non-relativistic "intense-field many-body S- 
matrix theory" (or IMST, review Q), in terms of the Dirac- Volkov-Feynman Green's function Gf{x,x') that satisfies: 

(ijt-e 4(u) - m)G F (x, x') = J 4 (x - x'). (10) 

A complete set of solution [?J of the Dirac- Volkov equation © can be conveniently written as: 

$(>)( X ) = . m*)ci + £^Hl)„w (ii) 

p V Po 2p- k p 



where, for a general pulse of transverse external electromagnetic wave, A(u), we get: 



5 ,eA(u) -p e 2 A 2 (u). , 
( o > du - 



(12) 



k ■ p 2k ■ p 

For a vector potential of a constant envelope it reduces to: f p (x) = —a p sin (u + \p) + b p sin (2m) + Q p u, with, a p — 

A °^cKrp ' = 8c^k~p cos d Xp — tan _1 [tan0 p tan(£/2)]; Up"-' are the four (two positive and two negative energy) Dirac 
plane wave spinors: 



^{2m{m + E) 
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HO = 



y/(2m(m + E) 



«;«,*= (3,4), (p„ = -£0 



(13) 
(14) 



where, E = + \/ {p 2 + rn 2 ) is always positive. The up and down positive energy 4-spinors tuW are given by: 



,(!=«) - 



and 



(2=d) 



Similarly, the up and down negative energy 4-spinors are: 



(3=u) 



and 



,(4=d) _ 
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Thus, the Dirac- Vokov Green's function satisfying the Feynman-Stuckelberg boundary condition (cf. Q), can be 
defined by: 



p s=l,2 



(17) 



s=3,4 



where, ^2 p — p~p / <^ 3 P- For an explicit momentum representation of G F (x, x'), which can be used conveniently for 
intense periodic electromagnetic fields, we refer to @, or Eq. (54) of @. 

Using G F (x,x') one may systematically solve Eq.® by iteration, to obtain the total Dirac state evolving from a 
specified initial state ip^ ( x ) '■ 



*,(a:) = i>i s) (x) + I d 4 x'G F (x,x')(e4(x'))^ s) (x') 



+ / dV I d A x'G F {x,x"){ lQ Vf{x")) 

J Xi J Xi 

xG F (x",x')(e4(x'))4 s \x') + 



(18) 



where Vf (x) is the external rest potential. This constitutes the Dirac wavefunction within the "relativistic intense-field 
many-body S-matrix theory" (or RIMST) subjected to the Feynman-Stuckelberg backward time positron interpre- 
tation of the negative energy "Dirac sea". It should be noted therefore that the present approach can be used for 
obtaining the probability of not only ionization, as we consider in this work, but also of such processes as pair creation 
in super-intense laser fields (cf. Q). 



III. GROUND STATE H+ MOLECULAR ION IN A DIRAC BASIS 



We shall represent the ground state four-component wavefunction ifj^ s \x), of hydrogenic molecular ions, in the spin 
state, ,s = (u,d), by a linear combination of relativistic atomic basis orbitals of a symmetry. To this end we define a 
four-component "Dirac basis" (with a symmetry) of the form: 

4 s) (r) = Bj(r)^(r>W;j = 1,2,3- •• . (19) 
R,(r) = N j T*- 1 e- x ' r (20) 

where rfjif) — 7 ■ n j{f), and the 4-vector ny(f) = (no,Tij(r)), with iiq — 1, and rijif) = iff^r, and /3j = m ^i rV .^ ', 
where, Vj and Xj are in general non-integer real parameters and Nj = (2Aj) t/j+1 ( 87rr 1 ^ : ' 2t/ ) 1 ^ 2 is a normalization 
constant. Note that the well known hydrogenic Dirac s-state wavefunctions Q in the present notation (from [lol ] p. 
418) can be reproduced by choosing a single Dirac basis state with V\ = Wl — (Za) 2 and Ai = mZa, where a is the 
fine structure constant and Z is the nuclear charge. The Fourier transform of the Dirac basis functions introduced 
above are evaluated analytically and we get: 



where, 



$\q) = / d\e-^^\r) = N jC °Aq)[l + g 3 {q)h • Q)W s) (21) 



and, 



- * T(1 + V lU %), (22) 
s°{q) = sin ((uj + 1) arctan (q/Xj)) , (23) 

9M = p'A - ^ v /r +(W ■ ]. (24) 

J q Vj V sin ({fj + 1) arctan {q/Xj)) 

The corresponding negative energy basis functions can be obtained from the F.T. (momentum representation) of the 
above positive energy basis functions by letting q — > —q, and replacing the positive energy up and down spinors 

, s = (1, 2), by the negative energy up and down spinors , s = (3, 4) defined above. Finally, we write the Dirac 

ground state wavefunction as: 

4 s \r;R)=^a j [Rj(r a ) V l j (f a )w^+Rj(r b )^(f b )w^} (25) 
j 

where, a,j are numerical constants, and the are the positive energy spinors defined above. 



IV. IONIZATION TRANSITION AMPLITUDE: UP AND DOWN SPIN CURRENTS 



Projecting the final Volkov state of 4-momentum p, Eq. (jlip . on the total wavefunction ^i(x), and retaining the 
leading significant term, we get the individual spin selected ionization amplitudes (s ~ u ot d to s' = u ot d): 

Ai on (s',s) = -i [ d i x^/\x)e4(x)4 s \x) (26) 



Using the Fourier transform of the ground state wavefunction, carrying out the d^x integration, and performing the 
somewhat lengthy Dirac 7-algebra, we arrive at the following explicit analytic result (in a.u. eo = mo = h = etc = 1): 



A lon (s -> s') = 2TTi^5(p -p b + (C P - n)K )T (n) (s', s) 



(27) 



with, 



2c 



(28) 
(29) 



Given the matrix elements M (s', s), the spin-specific rates of ionization become, 
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cpo\P\ 
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(30) 



The number of absorbed photons (denoted by n above, where hq is the threshold number) is determined by the 
4-momentum conservation relation, obtained from the delta function that appears naturally in Eq. (|28j) from the 



space-time integration: mo = e(, + ekin + Cp^i where ej = c c - y 1 c 2 — p\ ) is the binding energy (ionization potential) 



and ekin — c (^\/c 2 +p 2 — cj is the kinetic energy. On completing the Dirac algebra, the reduced Mj n \s', s)-matrix 
elements are found to be given by the following simple algebraic expressions: 



Mf'(u -» u) = B° n (mi + m29j{q)(p ■ q + i{p x q) z )) + B* ■ (m 2 p + m igj (q)q) 
-i (B* x (m 2 p - m 1 g J (q)q)) z , 



(31) 



M { f\u -yd) = m 2 g {q)B Q n * (i (p x q) x - (p x - i (B* X {m 2 p - m 1 g J (q)q)) i 
+ i B n x ( TO 2P - m igj {q)q)) v , 



(32) 



M [ f\d -> u) = m 2 gj(q)B°* (i (p x g) a |(px g) y ) - i (s* x (m 2 p - mi£,(g)g)) 
- x (™2P - m igj (q)q)) 



(33) 



Mj n) (d ->■ d) = £°* (mi + m 2 g J {q){p ■ q - i(p x q) z )) + B* n ■ {m 2 p + migj(q)q) 
+i (B* x (m 2 p - m igj (q)q)) z . 



(34) 



where, 



B° n = 
B„ = 

K ■ p = 
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(2J n + cosZ(J n+2 + J„- 2 )) 



Ack ■ p 

KO^O — K ■ p, K — (kq, K), Kq = Ol)/c . 



(35) 
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\Jn — \Jn {p>p i ^p; Xp) ^ ^ *7 n-\-2m {p>p)Jm {Pp)& 



i(n+2m)xp 



(38) 



are generalized Bessel functions of three arguments. The other parameters are the "field-dressed" electron momentum 
q = qq = p + (£ p - n) k, ( p = j^r^, e(£) = [e x cos(|) +ie 2 sin(§)], with (, = [0,±tt/2], p = V c 2 + p 2 = 

(n - C P ) k + y/c 2 -p 2 , m 1 = y/(p + c)/(2c), m 2 = -y/(po - c)/(2c), and iV po = 



Finally, the spin unresolved total ionization rate from an unpolarized target atom is obtained by simply adding the 
four spin-specific rates given above and dividing by 2 (for the average with respect to the two degenerate initial spin 
states): 
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(39) 



We may point out that Eqs. j28J), (|30| and (|39| that apply directly for the a g case, can be used to obtain also the 
probability of ionization from an analogous a u state, if one simply replaces the sum ("+") of the two interfering 
nuclear terms in these formulas with their difference ( "-" ) . 



A. Factorization of ionization probability and molecular two-slit interference 

It is interesting to see from Eqs. (|28j). (|30p and ([3"9"]) that the ionization amplitude and the probabilities factorize 
essentially into two factors, one of which depends on the nuclear geometry and the other on the electronic degrees of 
freedom. The geometrical factor leads to an interference effect arising from the superposition of the partial electronic 
waves that emerge from the two nuclear centers with a phase difference, q ■ (Rb — R a ) = q ■ R. This leads to the now 
well-known molecular two-slit interference (and its dependence on the symmetry of the electronic wavefunction of the 
active electron) first discussed in the non-relativistic case Q • Thus for example, for the ionization probability from the 
a g molecular orbital of Hj", the two-slit interference factor becomes: |1 + e lq ' R \ 2 which for q ■ R — 2mr, n = 0, 1, 2, ... 
gives rise to a constructive interference and an enhancement of the ionization probability by a maximum factor of 
4, while, for q ■ R = (2n + l)ir, n = 0,1, 2, a destructive interference yields a minimum in the energy-momentum 
distributions of the ionized electron. We note that the two-slit interference effect is generic in nature and can show 
its effect in all channels in which the active electron of the molecule determines the channel probability, not only 
in ionization but also e.g. in laser stimulated molecular electron diffraction, laser induced emission of (low or high) 
harmonic radiation etc. Furthermore, it can survive even in the signals that are obtained in orientation averaged 
measurements. In this case the interference probability factor reduces to: 

-l/^|l + e-T = 2(l + ^) (40) 

which indicates an interference modulation with the variation of the electron momentum q (or R), albeit decreasing 
in the modulation-depth with increasing momentum and/or the nuclear separation R. 



B. Linear and Circular Polarization, Photon Helicity and Spin Asymmetry 

The results derived above hold for any polarization of the laser field. We may note below the particular simplifica- 
tions that arise in the two most common cases of linear and circular polarization. 
1. Linear polarization: 

In the case of a linearly polarized laser field, the explicit formulas for the ionization currents, Eqs. (|30p and 
(|39|) . hold again when we simply put £ = 0, in them. Thus the only change occurs in replacing the generalized 
Bessel function of three arguments by the generalized Bessel function of two arguments: J n = J n (a p ,b p ), with 
o-p = a P {£ = 0),b p = 6 P (£ = 0). 



2. Circular polarization: 

Similarly, for the case of circular polarization of the laser field, the same expressions hold again when we simply put 
£ = ±7r/2 in them. In this case the generalized Bessel function of three arguments reduces to an ordinary Bessel 
function (of one argument) times a simple phase factor: J n = J„(a p )e ±m * p , with a p = a p (£ = ±7r/2), where ± stand 
for the positive and negative handedness of the circularly polarized laser field. 

Closer examinations of the expressions of the matrix elements M^ n \s' ,s) for circular polarization, as in the atomic 
case Q , show that the probability of spin-flip ionization for u — > d transition differ from the probability of the d — » u 
transition. It also suggests that the up and down spin-currents, even from a spin unpolarized (or equal mixture of up 
and down spin) ground state would show an asymmetry in any direction of emission of the electrons. This asymmetry 
is expected to persists even when the laser photons are of long- wavelength (or, the laser propagation vector k — > 0). 
This implies that the dominant component of the electron spin-current emitted from molecular targets by intense 
laser fields might be controlled by choosing the helicity of an incident circularly polarized laser field. 



V. SUMMARY 



In this paper, a 4-component Dirac relativistic theory of interaction of super-intense laser fields with molecular 
ion, and its heavier K-shcll analogs, is presented. To this end a "Dirac basis" of atomic orbitals is introduced and 
employed to construct the LCAOMO of the ground state of rlj. Explicit analytic expressions for the probability of 
ionization per second by an incident laser field of arbitrary intensity, wavelength and polarization is derived. The 
formulas derived apply for the total ionization probability as well as for the spin specific measurments of the ionization 
signal. They imply, among other things, an asymmetry of the probability of spin-flip, and of the up-spin and down- 
spin ionization currents, even from a spin unpolarized ground state molecular ion. The presence of a generic molecular 
two-slit interference effect, first noted in the non-relativistic case [H, is further confirmed in the present relativistic 
analysis. Finally, the possibility of controlling the dominant spin component of the electron current emitted from the 
inner-shell of a molecule by simply selecting the handedness of a circularly polarized incident laser field is pointed 
out. 
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